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Abstract 

It was argued recently that loop corrections to tree-level amplitudes are essential in the dis- 
cussion of the collisional energy loss of energetic partons in the hot quark gluon plasma: Instead 
of dE^ oll /dx ~ a 2 T 2 hi{ET/rn 2 D ), as first derived by Bjorken (assuming a constant a), the mean 
energy loss actually behaves as dE co \\/dx ~ a(m 2 D )T 2 . Here we calculate, within this resummation- 
improved framework, the probability distribution functions ('quenching weights') of the collisional 
energy loss. First results from a Monte Carlo implementation of this probabilistic collisional quench- 
ing shed new light on the interpretation of jet suppression in heavy ion collisions. 

PACS numbers: 12.38Mh, 25.75.-q 
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I. INTRODUCTION 



The dominance of radiative energy loss [1, 2, 3] for the phenomenon of jet quenching in 
heavy ion collisions seemed long established when Mustafa and Thoma [4] re-considered, in 
2003, the possibility of collisional partonic quenching [5, 6, 7]. However, only recently, after 
publication of new data on the nuclear suppression factor from non-photonic electrons [8] , 
this idea seems to have gained wider interest, see e. g. [9]. In fact, an additional contribution 
to the partonic energy loss appears necessary to reconcile data-adjusted model parameters 
with general expectations, say in the GLV formalism [3] the partonic rapidity density dN/ dy 
with the final entropy density, cf. [10]. Also the transport coefficients q inferred in [11, 12] 
within the BDMPS-Z formalism [1, 2] using probabilistic 'quenching weights' [13] are almost 
an order of magnitude larger than estimated [14]. In this situation, a collisional component 
of the energy loss is a welcomed remedy, as advocated in [9] . 

These (as well as previous) considerations rely on Bjorken-type formulae for the (mean) 
collisional energy loss, dE^ oll /dx ~ a 2 T 2 ln(£T / m 2 D ) . However, as pointed out recently 
[15], the collisional energy loss is an observable for which loop corrections to the tree level 
approximation are essential. Taking these higher order terms into consideration by using 
the running coupling in the calculation leads to modified predictions, which differ even 
parametrically from the commonly used expressions, see Eq. (3) below. 

In this note, in Section II, we will briefly outline the main idea of the (re-)calculation of 
the mean energy loss. In Section III, the probability distribution functions for the collisional 
energy loss are derived, which are necessary to describe jet quenching. These findings are 
implemented in a Monte Carlo model which will be presented, along with results, in Section 
IV. The implications of these exemplary studies with regard to heavy ion phenomenology 
are summarized in the Conclusions. 

II. MEAN COLLISIONAL ENERGY LOSS 

Following Bjorken [5], the mean collisional energy loss of an energetic parton (E ^> T) in 
a hot QGP (T ^> A) can be calculated from its weighted collision rate, which is determined 
by the flux and the cross section, 




(1) 
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The energy transfer in one collision, u = E — E', is related to the angle between the 'jet' 
j and the scatterer s as well as the invariant momentum transfer by t = —2(1 — cos8)ku. 
Furthermore, $ = 1 — cos#, and p s = d s n s is the density of the scattering partons, with 
d g = 16 and d q = 12n/ for gluons and quarks, and n±(k) = [exp(fc/T) ± l]" 1 in the ideal 
gas approximation. 

Since processes with small momentum exchange dominate, the cross sections can be 
approximated by the t-channel contribution, da/dt ~ a 2 /t 2 . In strictly tree-level approxi- 
mation, a is a constant parameter. Thus, the predictive power of the resulting expression is 
questionable because AE co n can probe a large range of the virtuality t where the running of 
the coupling, as a result of loop corrections, is important. On the other hand, these quantum 
corrections are also formally required to screen the long-range interactions in the medium. 

As recollected in [15], in a first approximation the loop corrections can be taken into 
account by evaluating the differential cross section with the running coupling 

where bo = (11 — |??,/)/(47r), and by imposing an IR cut-off in t, /i 2 ~ m|>, related to the 
Debye mass. Only the latter aspect has been utilized in the past, which led to Bjorken's 
formula [5] as well as various refinements thereof, which aimed at the precise determination 
of the cut-offs [6, 7]. However, as pointed out in [15], taking into account the running 
coupling modifies even the structure of the result. The t-integral in (1), with u ~ t, is 
basically J dt/{t\n 1 t) — — 1/lnt vs. Jdt/t — lntm the case of a constant coupling, hence 
the perturbative order and the dependence on the cut-offs is changed completely. The UV 
cut-off being related to the maximal possible energy transfer, tu max = O(E), the mean 
collisional energy loss becomes ^-independent for infinite parton energy [15], 

AE coll ~ a(/i 2 )T 2 Ax , (3) 

which is the direct imprint of asymptotic freedom. For a more detailed discussion of the 
relation of Eq. (3) to Bjorken's formula we refer to [16]. 

It was argued in [15] that such an 'improved' perturbative framework might be useful 
to estimate the energy loss even for rather small temperatures near T c as relevant in heavy 
ion phenomenology. The main idea is to realize that the energy loss probes a (complicated) 
QCD 4-point function in the special limit —t ^ s. The same 4-point function determines 



also (i) the distance dependence of the strong coupling as obtained from the heavy-quark 
potential at T = 0, and (ii) the screening of this potential at T > T c , as encoded in 
the Debye mass. Now both observables have been calculated in the large-coupling regime 
within lattice QCD. Remarkably, these non-perturbatively results - for both quantities - are 
reproduced quantitatively by the resummation-improved perturbation theory [17] with only 
one parameter, whose adjustment yields, for rij = 2, the canonical value Aqcd = 0.2GeV. 
This unexpected, but pleasing observation and the close relation of the three quantities may 
then support also an (otherwise not justifiable) extrapolation in case of the energy loss. 

In the revised form, the collisional energy loss is larger than previous estimates based on 
Bjorken-type formulae, e.g. AE q /Ax ~ 3GeV/fm for T = 3T C « 0.5GeV. It is emphasized 
again that, after fixing Aqcd from lattice QCD data, these results are basically parameter- 
free (while previous estimates had to 'guess' a value for the fixed parameter a.). 

III. COLLISION PROBABILITY DISTRIBUTION FUNCTION 

While the estimates above already hint at the importance of collisions for jet quenching, 
the mean energy loss is too crude to quantify the effect in heavy ion collisions. The fact that 
the initial hard parton spectrum is steeply decreasing requires to consider the probability 
density function [18], which in the context of the radiative energy loss is sometimes dubbed 
quenching weights [13]. 

For the collisional energy loss mechanism studied here, the probability density for an 
energy transfer uj to the medium along the path Ax reads 



Note that for larger Ax it is rather to be interpreted as the average number of collisions 



mean energy loss (1). 

Within the same calculation framework as outlined in the previous Section (for more 
details see [15]), the t-integral gives 



with the implicit constraint that |t*| = 2(1 — cos8)ku falls into the integration interval 





with an energy transfer uj. Obviously, its first moment 





da 

dt ti , 
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[/i 2 , 2(1 — cos 9)kuj max \. In the following we set cj max = E [23]; then P& x (w, T, E) vanishes at 
uj = E (and at higher energies). 

For the cross sections 'improved' by the running coupling in the form (2), the angular 
integral in Eq. (4) can be performed analytically and one readily arrives, without further 
approximations, at 

(5) 



for quarks and gluons, respectively. The dimensionless functions are defined by 

1(a) - I((b) 



Ja/b exp(C) ± 1 

with I(x) = xj \vl{x)— li(x), where li(rr) denotes the logarithmic integral function. In Eq. (5), 



X ± (a, b) are evaluated for a = /i /A and 6 = 4To>/A . In the following we will use 
[15], and determine the Debye mass self-consistently [17] 



2 - lm 2 
2 m D 



m 2 D 



(1 + |tv) 47ra(m!,)T 2 



(6) 



As depicted in Fig. I, P^(u;T) is, for typical temperatures, peaked at O(0.1)GeV. It is 
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FIG. 1: The probability density function (5) for collisions with an energy transfer uj, evaluated 
here for quarks, Ax = lfm and two representative temperatures. 

comforting to note that it decreases rapidly for larger values of uj, which guarantees that the 
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specific choice of cu max is unimportant, as anticipated above. For increasing temperatures, 
the distribution shifts to larger u and it becomes wider. 

For sufficiently small Ax, when P^ x {oj) is a probability distribution in the literal sense, 
the expression 1 — J duP^ x (uj;T) > is the chance that no collision happens within Ax. 
Accordingly, the zeroth moment of the probability distribution defines the mean free path 
A(T; E) of the hard parton in the medium, 

J duP x {u-T) = l. (7) 

It is instructive to estimate the size of the mean free length as a phenomenologically impor- 
tant quantity. First, one can verify again that predictions of the present approach are not 
very sensitive to the hard momentum scale; the mean free path lengths of a parton with 
intermediate energy differs from the asymptotic limit E — ► oo only moderately, see Fig. 2. 
Second, A is of the same order of magnitude as the particle distance, which can be estimated 




T (GeV) 

FIG. 2: The quark mean free path as a function of the temperature, for two different quark energies. 
Shown by the dotted line for comparison is the estimate of the typical particle distance tq in the 
medium. 

from the density p = ^j a d s n a - This is a clear sign of the strong coupling in the medium 
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(sQGP). While the quark free length coincides with ?"o = Po ~ T^ 1 at about T c , the 
gluon mean free path, due to the larger color charge, becomes smaller than Tq(T) already at 
T pa 0.5GeV. Interestingly, T ~ 4T C is close to the value estimated, by other means, in [19] 
for the 'transition' from a strong-coupling to a moderate-coupling regime. 

IV. PROBABILISTIC COLLISIONAL QUENCHING 

To quantify the effect of collisional energy loss, we develop a dynamical model of jet 
quenching in heavy ion collisions. Given the complexity of the phenomenon, which results 
in theoretical uncertainties on various levels, the intention is, for now, not to formulate the 
dynamics as realistically as possible; we rather focus on the transparency of the model. Hav- 
ing only few parameters, which hopefully capture essential aspects, might help to constrain 
the theoretical uncertainties in the underlying energy loss mechanisms. 

We consider here only the partonic energy loss due to collision within the dense plasma 
above T c , i.e., we explicitly 'switch off' radiation and we do not account for effects at and 
after hadronization. The first assumption is, of course, unrealistic but will provide for the 
benchmark we are aiming at. The second assumption, on the other hand, relies on the idea 
that a parton fragments only after it has left the medium. 

For the sake of simplicity, we study central collisions and focus on the mid-rapidity 
region. We assume a boost-invariant Bjorken expansion, T(t) = Tq(t /tq)^ 1 ^ 3 [2 I], with a 
homogeneous transverse profile. Consequently, the radius R of the longitudinally expanding 
cylinder is an effective parameter, say R pa 5fm [25]. The formation time of the QGP is 
related to the saturation scale, for RHIC To ~ 0.2fm is a generic value in the literature [26]. 
The initial temperature T , as the third parameter, is related to the total energy/entropy of 
the system; To ~ 0.5GeV is a representative value. To test the sensitivity of predictions on 
these parameters for the dynamics/geometry of the collision, they will be varied from the 
above values, as summarized in Table 1. 

The hard partons are initialized uniformly in the transverse plane (which adds to the fact 
that R is to be considered an effective radius). Depending on the 'jet vertex', characterized 
by r (< R) and the emission angle <p with respect to the radius vector, the partons propagate 
a length 

l± = (R 2 — r 2 sin 2 tp) ^ 2 ± \r cos 99) 



7 





set 1 


set 2 


set 3 


set 4 


set 5 


set 6 


R (fm) 


5 


3 


5 


5 


5 


7 


T (GeV) 


0.5 


0.5 


0.3 


0.7 


0.5 


0.5 


T (fm) 


0.2 


0.2 


0.2 


0.2 


0.1 


0.2 


t c (fm) 


6.1 


6.1 


1.3 


16.7 


3.1 


6.1 



TABLE I: The different sets of (effective) parameters to describe the collision. The last row gives 
t c = tq(Tq/T c ) 3 for T c = 0.16GeV, i. e. the QGP life time in the present approach. 



in the transverse plane (ignoring the small deflections by the collisions). By taking the 
modulus of the second term, l + > l_; below some results will be given separately for 'far 
side' and 'near side' jets. To segregate phenomena like the Cronin effect, we consider larger 
transverse momenta, pt <^ 4GeV, where the initial parton spectrum is calculable in pQCD; 
here we adopt the parameterization 

^{ 1 + 7J ' 

with Po w 1.75GeV and v w 8 [20]. 

In the course of time, as the hard partons propagate through the comoving medium, the 
Pt-spectrum evolves due to collisions with the thermal partons. The probability Pax(u; T) of 
an energy transfers uj along the stretch of path Ax depends on the local temperature of the 
rapidly cooling plasma. This important fact is taken into consideration in a Monte-Carlo 
simulation by discretizing the path such that Ax is smaller than the typical time scale on 
which Pax{w,T(t)) changes. Simultaneously, to adhere to the probability interpretation of 
Pax, Ax has to be smaller then the mean free length at the given time. For the implemen- 
tation, Ax = |A is an appropriate choice between numerical accuracy and computational 
efficiency 

Fig. 3 shows the result of a simulation, for set 1 of parameters, with 10 7 'jets'. The 
spectrum shifts towards smaller p t - although considerably different for the far and near 
side partons, as expected from comparing the life time, r c pa 6fm in this case, to the typical 
path lengths (obviously, /_ < R). These differences are better visible in the quenching ratio 



Qc 



dN^/dp t 
dN(T=o)/dp t 



(8) 
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FIG. 3: The evolution of the hard parton spectrum from formation time to freeze out, plotted 
separately for 'near-side' and 'far-side' partons, for set 1 of model parameters, cf. Table I. 

shown in Fig. 4. Although the present model is certainly too crude to make rigorous state- 
ments, such a distinct difference - even for central collisions and with no transverse expansion 
taken into account - is interesting with regard to recent results on jet correlations. Another 
noteworthy fact is the increase of Q q with p t - which is qualitatively different from previous 
predictions [4, 9]. The reason is that the later approaches were based on Bjorken-type for- 
mulae where, on average, alE^ oXX ~ \nE. However, as argued above, the dependence of the 
mean collisional loss on the parton energy is actually much milder: dE co \\ ~ E° for E — > oo. 

Finally, I briefly discuss the absolute size of the quenching factor. As discussed before, 
the collisional probability distribution function (5) is basically parameter free, so it is first 
comforting that, for a reasonable geometry/dynamics, the predictions without radiative loss 
do not undershoot the observed nuclear quenching factor, Raa ~ 0.2. Moreover, with an 
estimated light quark radiative quenching factor ~ 0.6 [9], a collisional quenching ~ 0.5 
(averaged over near and far side) seems just in the expected range [27]. Since the radiative 
modification factor is almost constant over a large pt-range [9], the moderate increase of 
Raa seen at large p t [8] could be a more indicative signature of the collisional contribution 
to the partonic energy loss. 

Before concluding it is necessary to investigate the sensitivity of the results on the model 
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Pt I GeV 

FIG. 4: The quenching factor corresponding to the spectrum in Fig. 3, given separately for 'near- 
side' and 'far-side'. For large pt events become rarer, which increases the statistical fluctuations. 

parameters describing the geometry/dynamics of the collision. From the results summarized 
in Fig. 5 it can be inferred that a very low initial temperature (set 3) leads to a diminished 
collisional quenching. A smaller effective radius (set 2) or a reduced life time due to a shorter 
formation time (set 5) give just a moderate change towards larger Q q values. A larger radius 
(set 6) only marginally decreases Q q compared to set 1, which is basically a consequence of 
the life time. On the other hand, larger initial temperatures give somewhat larger collisional 
quenching factors (set 4). 

V. CONCLUSIONS 

The fact that the average collisional energy loss per unit length has a different parametric 
behavior than widely held is reflected also in the corresponding probability density functions. 
In a resummation-improved framework, by taking into account the running coupling, these 
'collisional quenching weights' can be predicted unambiguously from the lattice-QCD ad- 
justable value of Aqcd ~ without having to assume a value for the coupling strength. Besides 
the mean energy loss, a second relevant moment of the probability distribution function is 
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4 8 12 16 20 

Pt I GeV 

FIG. 5: The quark quenching factor for the various sets of parameters summarized here by simple 
fits to the Monte Carlo results (depicted only for set 1) for pt 6 [4, 14] GeV (above 14GeV the 
statistical fluctuations become sizable). 

the mean free path of hard partons. Estimates for phenomenologically interesting tempera- 
tures yield the same order of magnitude, ^ lfm, as the typical particle distance in the QGP 
- which is another quantitative indication of a strongly coupled system. 

The probabilistic collisional energy loss can be readily implemented in a Monte Carlo 
code for jet quenching. A fully dynamical treatment of the quenching weights is feasible 
[28], thus a more involved dynamics of the rapidly cooling system can, in principle, be 
accurately taken into account. 

The present investigation is, in modeling the dynamics of heavy ion collisions, exemplary; 
more detailed results will be presented in a forthcoming study. Yet, because the main 
conclusion persists for a wide range of the principal dynamics/geometry parameters, it 
strongly suggests that the partonic energy loss has a sizable collisional component. 




set 1 
set 2 
set 3 
set 4 
set 5 
set 6 
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